Lifetime of metastable states in resonant tunneling structures 
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We investigate the transport of electrons through a double-barrier resonant-tunnehng structure 
in the regime where the current- voltage characteristics exhibit bistability. In this regime one of the 
states is metastable, and the system eventually switches from it to the stable state. We show that 
the mean switching time r grows exponentially as the voltage V across the device is tuned from the 
boundary value Vth into the bistable region. In samples of small area we find Inr oc \V — Vthf^^, 
while in larger samples Inr oc |V — Vth|. 

PACS numbers: 73.40.Gk, 73. 21. Ac, 73.50.Td 
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The problem of the decay of a metastable state has 
been addressed in a variety of areas including first-order 
phase transitions |l , Josephson junctions , field the- 
ory j^, magnetism ||], chemical kinetics Q. Meanwhile, 
progress in nanofabrication technology has made possi- 
ble observation of intrinsic bistabilities in double-barrier 
resonant-tunneling structures (DBRTS) ^ and superlat- 
tices . Recent experiments [g[ ^ with such devices have 
demonstrated that near the boundary of the bistable re- 
gion one of the two states is metastable, and its lifetime 
has been studied by measuring current as a function of 
time at different voltages. Thus, these devices provide 
an ideal experimental system for studying the decay of 
metastable states in real time. In this paper we develop 
the theory of switching times in double barrier structures, 
Fig. [^(a) . We expect the results to be relevant for other 
devices in which sequential resonant tunneling plays a 
key role in describing the electronic transport, such as 
weakly-coupled superlattices. 

We concentrate on the case of intrinsic bistability, 
which can be observed by measuring current / as a 
function of voltage V applied to the device while the 
impedance of the external circuit equals zero. As shown 
in Ref. g, for a certain range of bias V, two states of 





FIG. 1: (a) The potential profile of the DBRTS at applied 
bias V. (b) The I-V curve of the device has a bistable region 
between the dashed lines. The dotted line shows the process 
of switching from a metastable state to the stable one. 



current / are possible at the same value of the voltage, 
and the I-V curve has characteristic hysteretic behavior. 
As one increases bias, the upper branch ends at some 
boundary voltage Vth, shown schematically in Fig. |l|(b). 
If the voltage V is fixed just below the threshold Vth, the 
system stays in the upper state for a finite time r, before 
decaying to the stable lower state. 

We will show that the lifetime of the metastable state r 
can be understood by analogy to the problem of a Brow- 
nian particle in a double- well potential (Fig. ||). Here the 
coordinate of the Brownian particle has the meaning of 
the current / in the device (or the electron density n). In 
the problem of the Brownian particle, r depends expo- 
nentially on the height of the potential barrier Ub sepa- 
rating the local and global minima, i.e. r oc exp{Ub/T*), 
where T* is the temperature. Unlike a Brownian parti- 
cle, a DBRTS at nonzero bias is a non-equilibrium system 
in which fluctuation phenomena are driven by shot noise 
in the current rather than the electron temperature T. 
On the boundary of the bistable region, the local mini- 
mum disappears, and therefore Ub goes to zero. Thus, it 
is clear that r will depend exponentially on the voltage 
measured from the boundary Vth of the bistable region. 

Here we investigate effects of shot noise in DBRTS us- 
ing the framework of the theoretical model introduced in 
Ref. |lO[ The DBRTS is formed as a layered semiconduc- 
tor heterostructure. The electrostatic potential across 
the device is shown in Fig. |l|(a) . The potential is assumed 
to be independent of the x and y coordinates. The model 
includes only one subband in the quantum well. We fur- 
thermore assume that at zero bias the bottom of this 
subband Eq is above the Fermi energy Ep in the left 
and right leads. If the area of the sample S is small, 
we can assume that the charge in the well is distributed 
uniformly. Then, the state of the device is completely 
described by the electron density n in the quantum well. 
Below we will also discuss effects of non-uniform charge 
distribution in the well, which are important in the case 
of devices of large area. 

In the sequential tunneling approximation, the trans- 
port in the device is described by the following mas- 
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FIG. 2: Brownian particle in a double-well potential. The 
lifetime of the metastable state r depends exponentially on 
the height of the barrier Ub separating the local and global 
minima, r cx exp(f/i,/r*). The coordinate of the particle x 
has the meaning of the electron density n in the well, and 
U[n) = Su{n), see Eq. (||). 



ter equation for the time-dependent distribution function 
P{n,t) of the electron density n in the well, 
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Here /k, /p, and /q are the Fermi occupation numbers 
in the left lead, right lead, and the quantum well, re- 
spectively; Wkq(«) and Wqp{n) are the tunneling rates 
through the left and right barriers. The first two terms of 
Eq. (y) account for the processes which bring the system 
to the state of electron density n, and the last two terms 
describe the processes that take the system away from it. 
The first and the third terms on the right-hand side of 
Eq. describe tunneling of one electron into the well 
from the left lead, while the second and the fourth ones 
account for the probability of an electron in the quantum 
well to tunnel into the right lead. We dropped the terms 
describing the tunneling from the well to the left lead and 
tunneling from the right lead into the well. These con- 
tributions are negligible because the bistability emerges 
pof when the level in the well is close to the bottom of 
the conduction band in the left lead and far above the 
Fermi level in the right lead. (We assume T ^ Ep-) 

Assuming that the total number of particles in the well 
is large, nS 3> 1, we can expand in the small parame- 
ter 1/ S. Keeping terms up to the second order we reduce 
the master equation to the Fokker-Planck equation. 
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The exact expressions for A{n) and B{n) are rather com- 
plicated, but near the threshold they can be calculated 
analytically, see Eq. below. 

The stationary solution of Eq. can be easily ob- 
tained: 
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where u{n) is the effective potential. 

In the derivation of the Fokker-Planck equation (||), 
the coefficients A{n) and B{n) appeared as the first and 
second terms of the expansion in 1/ S. We therefore con- 
clude that A{n) / B{ti) cx S, and is independent of S. 
Thus when the area S is large, the distribution function 
Po("-) has very narrow peaks near the minima of u{n). 
If we neglect the fact that the width is finite, then the 
electron density n in the well can be found by minimizing 
u{n). The minimization condition written as A{n) = is 
in agreement with the results of Ref. 

Our model allows for an analytical treatment at small 
values of the parameter A = me^/27r?i^C, where C is 
the capacitance of the device per unit area. Then the 
calculations are greatly simplified, and one obtains the 
following expressions for A and B near the threshold: 
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Here 5V = Vth — V, the electron density at the threshold 



nth 



^ ( ^ ) pf ; and T^ ji are the transmission co- 



efficients of the left and right barriers at energy Eq and 
zero applied bias. If A is not small, we cannot get ex- 
plicit expressions for a, 7 and 6, but the generic form of 
Eq. (^ ) remains unchanged. 

The potential u{n) is shown schematically in Fig. || for 
a voltage which lies slightly below the threshold voltage 
Vth- Close to the threshold the potential can be approx- 
imated by a cubic polynomial, 

u{n) K, -a{n - nth) + ~ ^t^ f + u{nth)- (5) 

Taking into account the exponential dependence of the 
mean switching time on the barrier height, r — tqc^'', 
and using Eq. (ph, we find 
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(6) 



The prefactor tq can be found using the techniques de- 
scribed, e.g., in Ref. 
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It is important to note that the form (g) of the poten- 
tial u{n) and the hnear dependence a oc 5V are dictated 
by analyticity of the potential near the threshold. Thus, 
the applicability of the following results is not limited to 
a particular model of transport in DBRTS. A 3/2-power 
law analogous to Eq. ^ was theoretically predicted for 
different physical systems in Refs. ||^, ^, ||, Experi- 
mentally it was observed recently for an optically trapped 
Brownian particle [|l3j . 

The result ^ has been obtained under the assumption 
that the electrons spread rapidly in the x-y plane, and 
their density n is uniform. In large samples, however, the 
spreading takes a long time, and one has to account for 
the dependence of the density n on the point r = (x, y) 
in the well. This can be done by generalizing the Fokker- 
Planck equation (||) to the case of non-uniform n{r). 

We begin by studying the in-planc diffusion of electrons 
in the well neglecting coupling to the leads. For simplic- 
ity we neglect the electron correlation effects; the inter- 
actions of electrons will be accounted for in the charging 
energy approximation. Assuming that the electrons dif- 
fuse independently, one can write a master equation for 
the distribution function P{n{r),t} as follows. During 
the time At at most one electron can move from position 
ri to r2, that is. 



P{7i(r), t + At}- P{n{r),t} = J J dndra 

X P{n{r) + <5(r - n) ~ S{r - r2),t}W{ri,r2, At) 



-P{n{r),t}W{r2,Vi,At) 



(7) 



Here M^(ri,r2, At) is the probability density of an elec- 
tron diffusing from a point ri in the plane of the quantum 
well to point r2 during the time interval At. Since elec- 
trons are fermions, a particle can diffuse only from a filled 
state at ri to an empty state at r2. Assuming that the 
diffusion is due to the elastic scattering of electrons by 
defects, we find 

W(ri,r2,At)-5(ri-r2,At)j. J ME)[1 - h{E)] dE, 

(8) 

where u is the density of states in the well (per unit area), 
and fi.2{E) are the occupation numbers at points ri, r2. 
The classical diffusion probability g(r, At) is given by 
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where D is the diffusion coefficient. The approximate 
form is obtained in the limit At 0. 

Using Eqs. (^), (||) and expanding the distribution 
function from Eq. up to the second order in 5{r—ri) — 



5{v — Y2) we obtain a functional Fokker-Planck equation, 
dP{n{v),t} 
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= -(A'i-A'2)/T 



Here /ii and /i2 are the electrochemical potentials at ri 
and r2 , respectively. Their values are found by adding the 
electrostatic potential e^n/C to the Fermi energy n/z/. 



Mi,2 = -~-n(ri,2). 



(11) 



Here C is defined by e^/C = e^/C + 1/v. 

Substituting Eq. (|ll|) into Eq. (|l^, integrating twice 
by parts and assuming that |^i — /i2| <C T, one obtains 
the following Fokker-Planck equation. 
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(12) 

The stationary solution of Eq. ( p^ is found easily by re- 
quiring the part of the integrand after the first functional 
derivative to vanish. 



Po{n} = exp 



e^n^(r) 
2C 



dr 



This result has a simple physical meaning. The equilib- 
rium distribution function Pq has the Gibbs form e"^/"^, 
with the energy per unit area e^n^/2C in agreement with 
the electrochemical potential (|ll|). 

Using Eq. ( p^ ) we can take into account processes of 
charge spreading in the quantum well. The electron den- 
sity n(r) can change due to either tunneling of electrons 
through the barriers or their diffusion inside the well. 
Thus, we must add the terms from the right-hand side of 
Eq. (||) to Eq. ( p^ to account for both processes. The 
combined Fokker-Planck equation takes the form 
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Here we neglected the second term in Eq. ( p^ . This can 
be done as long as the temperature of the electrons in 
the well is much lower than the Fermi energy [T5| . 
The stationary solution of this equation is 
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const 
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(14) 
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y^ajj Zs{r/ro)- Therefore, one can calculate the loga- 
rithm of the mean switching time r as F{ns} — F{nmin}-, 



FIG. 3: Numerical solution of Eq. 

where t] = vDe^ /CSB{n). Note that in the limit D oo 
the electron density is uniform, Vn = 0, and we recover 
the result (1). 

In the case when the diffusion coefficient D is finite, 
the electron density varies from point to point in the 
quantum well; hence, this problem is infinite dimensional. 
In the multi-dimensional case, the system escapes from 
the local minimum of potential through a point where the 
barrier separating it from the global minimum takes the 
lowest possible value, i.e., through a saddle point. The 
mean switching time r is determined by the potential 
at the saddle point, measured from the local minimum. 
This approach is similar to the one used in the theory of 
kinetics of first-order phase transitions , with F playing 
the role of the free energy. 

The saddle point of the functional F{n} is achieved at 
n = Us (r) which almost everywhere in the sample is very 
close to the density rimin of the system at the local min- 
imum of F . However, in a region of some characteristic 
size ro, the density nsfr) changes in the direction of the 
global minimum. Fig. g. Thus, the DBRTS of large area 
first switches to the stable state in a region of size rg, 
which then expands to the whole sample. 

We perform the following calculations in the regime 
of voltages very close to the threshold Vth, where we 
can use u(n) in the form (|5|), and rj takes a constant 
value rj = vDe^ /Ch. Initially the system is in the lo- 
cal minimum, described by a uniform density rimin = 
iT'th + \/ ct/l- In order to find the saddle point, it 
is convenient to parametrize the electron density n(r) 
in terms of a dimcnsionless function z{p)^ such that 
n(r) = rimin — yfapf z{r / tq) , where the length scale 
ro = (a7)~^/'*. Using this parametrization, one can 
find the saddle point of F as a non-trivial solution Zs{p) 
of the equation 



0, 

: at p - 



(15) 
It can 



with the boundary condition Zg (p) 
be obtained numerically, Fig. |[ 

The switching rate is proportional to the dis- 
tribution function Po{ns}, where ns{r) — nmin — 



1 ^'^ AT/ 
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Here the constant k 
According to Eq. 



62.01 was found numerically |16 



If 



In ^ does not depend on the 



area S. On the other hand, since the critical fluctuation 
ns(r) can be centered anywhere in the sample, the switch- 
ing rate is proportional to the area of the sample S*, 
hence T2 csc 1/5'. The exact calculation of the prefactor 
T2 presents a number of theoretical challenges, which we 
leave for future work. 

In contrast to the case of small samples Eq. (^), the 
logarithm of the escape time (|l6|) in large samples is 
linear in 6V . The crossover between the two regimes 
occurs when the area S of the sample is of the order 
277(q!7)^^/^. One can see from Eq. (^) that 
5^-1/2. Thus, one can observe this crossover in a 
single sample by tuning the voltage. Indeed, at relatively 
small 5V we will have S* ^ r§ and Inr oc SV^/"^, Eq. (13), 



of rl 



whereas at larger 5V we have Inr oc 5V , Eq. (|16[). 

In conclusion, we have studied the switching time t 
from the metastable state to the stable one in DBRTS. 
We showed that r is exponential in the voltage measured 
from the boundary of the bistable region; it is given by 
Eq. (H) or (^6|) depending on the area of the sample. Our 
results can be tested in experiments similar to Refs. |[ |^. 
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